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ABSTRACT 


Transient film condensation on horizontal surfaces 
and laminar film condensation flow on a horizontal isothermal 
surface are considered. Transient film condensation on hori- 
zontal surfaces is analyzed, for various boundary conditions 
of the plate. These ares (a) isothermal surface, (b) plate 
with insulated, bottom surface and a negligible temperature 
variation across its thickness, (c) semi-infinite approxi- 
mations, (d) plate of finite thickness and bottom surface 
at a fixed temperature, and (e) plate of finite thickness 
and bottom surface insulated. Each case is considered for 
two inviscid flow conditions, one, a specified fraction of 
the rate of condensation, occurring at any time, flowing out 
across the edges of the condensing surface at that instant, 
and, second, the outflow d omina ted by gravity. The appropriate 
conservation equations are simplified by using several approxi- 
mations valid for many physical circumstances of interest, and 
the resulting ordinary non-linear differential equations, 
except for the first type of flow in case (b), whore exact 
solution is obtained, are integrated numerically. The theore- 
tical results are shown to depend on two dimensionless para- 
meters, one for the condensation, except for the case of 
isothermal surface, and, the other for the flow. The situation 
of no-flow arises when the flow parameter is zero . 



Laminar film condensation flow on * horizontal 
isothermal surface is analyzed in three singes (a) transient 
film condensation without flow, (h) transient laminar flow, 
and (c) steady, laminar flow considering the surface tension 
effects, fourth order non-linear differential equations are 
obtained for transient and steady flows, negligible surface 
tension, or large values of Bond lumber reduce them to 
second order. The condensate film profiles for steady laminar 
flows are found to be strong functions of edge conditions at 
the end and the length of the condensing surface, and arc 
also observed to depend on the vapour properties and the 
difference between the temperature of the vapour and that of 
the condensing surface. Results obtained show that for the 
condensing surface, terminated by a round foil, the surface 
tension effects are negligible, even for very small values 
of the Bond number. The analysis and results for transient, 
laminar flow show that the transient effects are large for 
highly viscous condensing fluids, or for low gravity field. 

A very, simple expression for the local Hus sol t number is 
also obtained. 
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CHAPTER 


I 


Tlie importance of the condensation heat transfer process 
is due to the large number of applications it finds in industry* 
Dor example, the condensation process is of considerable 
importance in the design and study of the condensers of 
thermal power stations, thermal processing systems, nuclear 
power generation systems, refrigeration and air-conditioning 
systems and in many manufacturing processes, which employ the 
condensation of a hot vapour for the thermal treatment process. 
Condensation soldering, curing, heat treatment and degreasing 
processes are few examples of the applic tion of condensation 
in manufacturing industry. 

A vapour in contact with a surface at a temperature 
lower than the saturation vapour temperature, corresponding 
to the given vapour pressure, will condense on the surface. 

The condensate thus formed will be in the subcooled state due 
to contact with the cooled surface, and • more vapour will con- 
dense on the exposed surface and on the previously formed 
condensate. Depending on the behaviour of the condensate, with 
respect to the cooled surface, the condensation process has 
been zdivided into two distinct condensation mechanisms. If the 
condensate tends to wet the surface, thereby forming a liquid 




film, the process is termed film condensation. If the conden- 
sate does not wet the surface, out, instead, collects in growing 
droplets on the cooled surface, the process is termed dropwise 
condensation. Both types of condensation are common, and, in 
some circumstances, may be simultaneously encountered on a 
single surface. Due to the direct contact between the vapour 
and the surface and the extended - surface effect of the 
droplets the heat transfer rates are much larger in the case 
of dropwise condensation, as compared to those in film conden- 
sation, where the film introduces a thermal resistance . 

If the condensate forms a film over the cooling surface, 
the process is quite different, file heat transferred to the 
cooled surface, from the vapour, occurs by the convective heat 
transfer across the liquid film, the higher temperature being 
at the interface, which is in contact with the vapour. ‘Phis 
convection process may greatly reduce the overall heat transfer 
rate, depending upon the tb.iclcness of the film, the fluid 
properties and turbulence in the flow. Therefore, the heat 
transfer rates for film condensation arc considerably lower 
than those for dropwise condensation, and in many -applications, 
efforts are made to generate a dropwise condensation by 
vibrating the cooled surface. 

Host of the condensation processes, generally encountered, 
give rise to film condensation or are mixed in character. Since 
a developed, steady, laminar film will also give the lowest 



steady lieat flux, one would have to contend with in any 
design-, a study of the film- condensation process takes on a 
special importance. Almost an;* departure from the simple 
film condensation process will give rise to a higher heat flux, 
at a given temperature difference between the surface and the 
vapour, At. In transient processes the heat transfer 
coefficient varies with time and may give a minimum, value as 
observed in transient natural convection flows. 

The physical nature of the s testy laminar film condensa- 
tion process is well understood, duo to the many previous and 
recent investigations. The steady laminar process was first 
considered in 1916 by hum soli [l] and his predictions were 


found to be in good agreement with experimental results for 
condensing fluids of Prandtl Humber, lip? around unity. This 
simple analysis ignored acceleration, or inertia, effects, 
assumed a linear temperature distribution in the film and 
neglected energy effects with liquid subcooling and interfacial 
shear between the liquid film and the vapour. Subsequent 


analysis by Brcomley [2], 


and by fobs enow [ 3 ] dispensed with 


some of these approximations. The results are identical with 
those of the Ifusselt theory , except that h- in those relations 

X £5 


is replaced by h^ + g- A t; where C 
specific heat of the condensate liquid 


P 


and h-> a. re the 
1 l!i 

and the latent heat of 


the vapour, respectively. The correction 
negligible for some condensation processes 


C. At is not 

o 

and gives results 
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which, deviate from those obtained by Nusselt [l] for large 

values of C At/ h„ r .. In practice, this value is usually 

Pp - I -'i 

less than 0.1 [4]. 

A boundary la^/er similarity analysis, neglecting the 
interfacial shear, was carried out by Sparrow and Gregg [4,5] > 
without making anjr of the other three as sumptions of Musselt, 
listed above. Por fraud tl Humber unity and above and 
C ht/h ± 1, the results obtained by Sparrow and Gregg [4,5] 

P n I 

are quite close to those obtained by omitting the inertia terms 
But inertia terms play a significant role in the condensation 
heat transfer process for Prandtl Humber values much less than 
unity, generally for Pr^ < 0.05. The results of ITusselt [1], 
Rolls enow [3] and Sparrow and Gregg [4] have been compared 
graphically by Gebhart [6]. 


Other film condensation analyses, considering the 
shear at the liquid - vapour interface, were presented by 
Koh et al.[7], Chen [8] and Poors and Ililes [9]. The results 
show that the effect is very small for Prandtl lumbers of 10 


or greater. The effect increases with the subcooling parame- 
ter 0^ A t/h^ and, for very low Pr-^, the reduction in heat 
transfer may be as large as 50 percent, for C At/h. 0fv as 
large as 0.1. Rohsenow and Hartnett [25a] have presented the 
comparison of results of Chen [8] and Sparrow and Gregg [4,5]. 
Yang [10] used a series technique for film condensation with 


variable surface temperature and obtained solutions for a 
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range of Prandtl Humber. Sparrow and Ian [ll] analyzed the 
effects of noncondensable gases present in the condensable 
vapour and reported that the presence of a very small amount 
of noncondensablo gas reduces the heat transfer coefficient 
significantly. Ilinkowycz and Sparrow [12] considered the 
effect of superheated vapour and have reported that the rate 
of condensation is slightly increased by the superheat in 
vapour. Sukhatme and Rons enow [13] and Mills and Sehan [14] 
have shown that the interfacial resistance at the vapour 
liquid interface is not important for many fluids and practical 
conditions. Steady, laminar film condensation, in the 
absence of the gravitational field, has also been investigated 
by 3Ioh [15], Sparrow et al. [16], Minkowycz and Sparrow [17] 
and G-ebhart [ 19 ] . 

There is apparently much less in literature on the 
transient flow and on laminar film condensation arising on 
horizontal surfaces. These situations are typically much more 
complicated than those described above. However, transient 
and steady film condensation on finite horizontal surfaces are 
important in technology and in nature, as, for instance, during 
the start-up of systems and in a number of other processes of 
interest where the horizontal surface is immersed in a conden- 
sing vapour region at a given time. 

Sparrow and Siegel [13] were the first to analyze the 
transient film condensation flow on vertical isotherma.1 
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surfaces. Their analysis ignores inertia effects and assumes 
a linear temperature profile in the condensate film. Himmo 
and leppert [20, 21] have analjrzed the laminar film conden- 
sation flow over horizontal isothermal surfaces neglecting the 
inertia, momentum and surface tension effects, and assuming 
a linear temperature distribution in the condensate film and 
zero shear at the vapour - liquid interface. 

Gebhart [19] and Rohsennow and Hartnett [23b] have 
considered the transient film condensation over isothermal 
horizontal surfaces with no run-off, i.e., no flow of con- 
densate over the edges, and have found that the condensate 

l/2 

film thickness 5 is given by, 6 = iC-j % , where t is the 

time and a constant, which depends on the vapour and 
condensate liquid properties and on At. The latter investiga- 
tors have also analyzed the transient condensation problem, 
with zero gravity. All these analyses are for a uniform trail 
temperature, assuming a saturated va.pour and neglecting the 
resistance at the liquid -va.pour interfa.ee, the momentum a.nd 
convection effects and the variation of fluid properties with 
tempera.ture. 

Pfahl [22] has shorni that some of the conduction 
solutions, available in Cars law and Jaeger [26], for problems 
with a change of phase, can be employed in the modelling of 
liquid - vapour condensation problems, provided that there is 
no motion in the liquid phase. In particular, 


he has 
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considered two exact solutions, for transient film condensation 

on horizontal surfaces ; one, on an isothermal surface, and 

the other, on a semi-infinite container bottom, initially, at 

temperature t^ < t g . In both the cases, the condensate film 

l/2 

thickness 5 increases with time and is given by, 6 = 1 2 t , 
where K 2 is a constant and depends on the vapour, condensate 
and surface properties and on At. Recently, Mollendorf and 
Chu [24] have analyzed the transient film •condensation over a 
horizontal body, initially at t- (< t_), for an insulated 
underside and a low value of the Biot Slumber, i.e., negligible 
temperature variation ::across the thickness of the body. In 
their analysis, they have considered the gravity ~ dominated 
flow of the condensate. Their analysis is also subject to 
the usual approximations, mentioned above, and neglects the 
surface tension effects. 

The present work is directed at an analysis of film 
condensation on horizontal surfaces. The work is divided 
into two main parts. The first concerns the transient film 
condensation on a horizontal body of uniform thickness, with 
run-off, and the second concerns the laminar film condensation 
flow on a horizontal isothermal surface. 

For the transient film condensation on horizontal 
aolates, the analysis has been carried out for five different 
cases. These are related to flow and heat transfer for the 
following conditions, assuming condensation to occur on the 
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upper surface: 

(i) An isothermal surface 

(ii) A plate with insulated bottom surface and a 
negligible temperature variation across its 
thickness. 

(iii) Semi-infinite Approximation: Solution for 
short times and thick plates. 

(iv) A plate with a finite thickness and bottom 
surface maintained at a constant temperature. 

(v) A plate with a finite thickness and an insulated 
lower surface. 

Furthermore, each case has been analyzed for two types of flow 

I. A specified fraction of the rate of condensation, 
occurring at any time, assumed to flow out across the edges 
of the condensing surface. 

II. Gravity dominated flow, in which the outflow velocity 
is determined by the gravitational force resulting from the 
height of the condensed fluid. 

The results obtained show that the heat transfer 
coefficient is a strong function of time and is infinite 
at time x = 0, and decreases exponentially to zero for large 
t. A dimensionless condensation parameter (3, depending on 
the fluid and the plate properties, (t fj - t ^ ) and h^ , has 
been found, which controls film ■ condensation on horizontal 


A 



9 


plates, for all boundary conditions considered except for the 
isothermal surface. For no flow and for a specified fraction 
of the condensate flowing out, the condensate film thickness 


increases with time except for the cases in "hich lower surfa.ee 
of the plate is insulated. Por the plate with insulated 
bottom surface, the film thickness achieves a constant value 
after a large time, when the temperature of the plate has 
reached a value close to the saturation temperature of the 


vapour. In the case of the semi-infinite solution, the 
temperature of the upper surface of the plate picks up a 
constant value, depending on the fluid and the plate properties 
as soon as it is exposed to the condensing vapour. In all 
other cases the temperature in the plate is a. strong function 
of time except for the isothermal surface, 'For no flow condi- 
tion, the results obtained for an isothermal surface and for 
thick plates are in good agreement with the exact solutions 
presented by Pfahl [22]. 


Por gravity dominated flow, another dimensionless flow 
parameter y has been obtained, which, along with ;3, describes 
the condensation process over horizontal plates. In the 
present case, the condensate film thickness first increases 
with time and then becomes constant, the rate of condensation 
being equal to the rate of flow, except for the cases of 
insulated lower surface, where the film thickness becomes 
almost zero, the surface temperature being very close to t,,. 



10 


In the second part of the present work, an analysis 
of the laminar film condensation flow on a horizontal 
isothermal surface has been carried out. The following three 
stages have been considered. 

1. Transient film condensation without flow. 

2. Transient, laminar film ■ condensation flow. 

3. Steady, laminar film condensation flow. 

Non-linear differential equations of fourth order have 

been obtained for transient and steady state laminar film 
condensation flow with appreciable surface tension effects. 
For large values of the Bond Number, Bo, or negligible 
surface tension, they reduce to second order systems. The 
differential equation obtained for steady laminar flow with 
large Bo is the same as the equation obtained by Hiramo and 
Leppert [20, 2l], 

For st eady , laminar flow, the results obtained show 
that the condensate film profile is a strong function of 
edge condition, i.e., P min > the dimensionless film thickness 
at the end, and the length of the condensing surface ,1'. khile 
the rate of condensation, and, hence, the flow rate increases 
as L is increased, they decrease with increased r) min . Tke 
parameter K, which is a function of the properties of the 
condensing va.pour and At, also affects the film profile. The 
condensate film is found to be thicker for a higher value of 
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K, and, hence, reduces the flow rate. The local Fuss cl t 
Number Nu is equal to x/6(x) where x is the location and 

X. 

5(x) is the film thickness at x [Pig. 22], 


Por transient laminar flow, a parameter R has been 
obtained which is largo for highly viscous condensing 
fluids, or for low gravity, and is found to increase with 
the length of the condensing surface, L, and C At/ hi . 
The rate of increase in film thickness decreases with time 
and, hence, the local heat transfer coefficient h__ is also 

-X 

reduced. The surface tension effects are found to be 
negligible for the condensing surface terminated by a round 


fall. 


All the computations are performed on IBM 7044 computer 
at I IT Kanpur. Por the numerical solution of ordinary 
differential equations, for transient film condensation, 
combined with conduction, on a horizontal body (Chapter II), 
4th order Runge-Kut ta method is used and for the ordinary 
and partial differential equations, developed for the laminar 
film condensation flow, the finite - difference technique 
is employed. 



CHAPTER II 


TRAITS IE NT FUJI COMDSITSATIOH, COISIHED RITE COITDUCTIOIT, 

OH A HOPJiOPIAI. BODY 


2.1 Analys is ' 


If the upper surface of a horizontal body, initially 
at temperature t. , is subjected to a saturated vapour, at a 
uniform temperature t 0 , corresponding to its pressure, 
condensation will take place on the surface provided t^ < t . 
The rate of condensation depends on the imposed boundary 
conditions for the body, besides other factors , such as fluid 


properties, latent heat of vapour, temperatures etc. This 
chapter considers the analysis and the results for transient 
film condensation on a horizontal surface with flow across 
its edges. The complete analysis has been divided into five 
parts, each part related to a particular circumstance of the 
imposed boundary conditions, for the horizontal bocljr of 
uniform thickness. All these different cases are often of 
interest in many practical problems related to 'transient film- 
condensation on horizontal surfaces. These cases have already 
been mentioned in the previous chapter. 


The analysis proceeds by applying the principle of 


mass continuity 
of energy to the 


to the condensate .film; and that of conservation 
horizontal plate. A similar analysis for the 



steady state film condensation on a vertical isothermal 
surface was first carried out by Husselt [1]. Several of 
the results obtained from the boundary layer analysis of 
Sparrow and Gregg [ 4 ] » for steady, laminar film condensation 
on a vertical isothermal surface, are of interest to the 
present work. In particular, the analysis presented in this 
chapter is subject to the following approximations, besides 
any other mentioned in the course of analysis. 

(a) Ho interface resistance at vapour-condensate interface, 
i.e. , the vapour-liquid interface is at the vapour saturation 
temperature . 

(b) The interfacial shear between the liquid film and the 
vapour is negligible. 

(c) The vapour is saturated and is at a constant tempera- 
ture corresponding to the pressure. 

(d ) The physical properties of the liquid and the vapour 
are constant with respect to temperature variation. 

(e) Ho surface resistance to heat transfer at the 
c ond e ns ing surf ac e . 

(f) Negligible inertia and surface tension effects. 

(g) The temperature distribution in the condensate film 
is linear. 

The assumptions (a) to (f) are ver;?- common in conden- 
sation analysis [1-5]. Regarding the last assumption, it 



lias been shown by Rohsenow [3] that, for the range of 
parameters encountered in most applications, a linear tempera- 
ture profile serves as an excellent approximation. Also, 
Sparrow and Gregg [ 4 ] have shown that the film temperature 


distributions are essentiall' 
CL At/li , in steady flows 
findings can be carried over 
be inferred from the results 
t ra ns i ent i c e -f o rma t ion [30] 


■ linear for email values of 
si hat all these steady state 
to the transient situation can 
.from a study of the problem of 


2.1.1 Isothermal Surface : 


ibis is a very simple situation, of transient film 
condensation, in which a surface maintained at 0 . temperature 
t^ is exposed to a saturated vapour at a constant temperature 
t_ at time t = 0 (Pig. 1). H ith the approximations a. .so; 
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fig. It Schematic showing the physical model of 

transient film •condensation on a horizontal 
is o thermal surface . 


oar-lior 


XX1C 


problem is analyzed for the two 


.disoussod 
types of flow. 

first, a specified fraction of the rate of condensation 
at any time, flowing out across the edges of the 


occurring 
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condensing surface is considered. It is assumed that the 
condensate which flows out is at the saturation temperature 
t g . The mass of vapour being condensed cat any time is, thus > 
given by: 



where 6 is 
A xlie area 


the condensate film thickness at 
of the condensing surface, P the 


any gi ve n t in e , 
density of the 


condensate liquid, and r the 
condensed per unit time, at a 
the edges of the condensing s 


fraction of the condensate, 
given time, flowing out across 
urfa.ce of the body at that time. 


Assuming that the latent heat released at the vapour- 
liquid interface is conducted through the condensate film 
without any loss, we obtain, 


P h. 

Ti 


& . pA L x 1 

r J dt 


v 

1 


t 


~&(%T 


X 


(2.2) 


Integrating the above 
at time x = 0, 6=0, 


equation and using the condition that 
since condensation starts, at x = 0, we 


get. 


O 




M 


2h-j (t r . - t.) 

— ' L - , - J - 1 - - (1-r) 

P ll.,.. 


(2.5) 


for the no flow condition, i.e., r = 0, the Bq. 
the same as that obtained by G-ebhart [15], and 


(2.5) is 
by .Rolls enow 


and Hartnett [25b]. 
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condensing surface is considered. It is assumed that the 
condensate which, flows out is at the saturation temperature 
t g . The mass of vapour being condensed at any time is, thus 
given by: 



where 6 is 
A the area 


condensate 


the condensate film thickness at any given time, 
of the condensing surface, P the density of the 
liquid, and r the fraction of the condensate , 


condensed per unit time, at a given 
the edges of the condensing surface 

Assuming that the latent heat 
liquid interface is conducted throng 
without any loss, we obtain, 


time , flowing out across 
of the body at that time 

released at the vapour- 
h the condensate film 


P h.„„ - c / \ 

x g o o ( x ) 

Tl~r) 'dr 


— V 


x i £ - J h'(xT' ] 


( 2 . 2 ) 


Integrating the above equation and using the condition that 
at time t = 0, 6 


+ - Sr ' 1 - ^ = o . n = o, since condensation starts, at % = 0, we 


get : 


S(t) = 


M 


SJV-.V. (1 . r) 

P h 

h ± V 


(2.5) 


for the no flow condition, i.e., r = 0, the 3q. (2.5) is 
the same as that obtained by Gebhart [19], and by Rohs enow 
and Hartnett [23b], 



we define the 


Now, for generalization of the anal^rsis, 
temperature excess 9 as; 

0 = t - t^ (2.4) 

and non-cl inensionalize 5 and t with the characteristic length 
6 C and time v , respectively, as; 

6 


Dimensionless film thichness, \ = 

and Dimensionless lime, E = ---• 

T c 

Here , 


(2.5) 

( 2 . 6 ) 


2 v ft 
bf ^rH3_ 

\j a phj, 

xg 

o 


a 


(2.7) 

(2.3) 


where b and a arc the thickness and the thermal diffusivity 
of the plate. 

It may be noted that for an isothermal surface, the 
variables b and a do not enter the analysis. However, this 
non-dimensionalization has certain advantages in the other 
cases to be analyzed later. In the present case, i.e. , 
for an isothermal surface, we can use an arbitral value for 
b 2 /a, 1.0 being the simplest. But if a comparison of the 
results of this case is to be made with those for other cases 
to be discussed later, it would be convenient to take the 

o • 

anno. Hith the above definitions , we may now 


2 / 

am value for b / C H 

so thi 

ales arc the sane. 

Hith 

■it o Eq. . (2.3) as 
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x = ^2(1 - r) - 5 


(2.9) 


It is also important to consider the transient film 
condensation on an isothermal surface, for the second case 
when the flow is gravity dominated. A balance between the 
heat released in the condensation process and the energy 
conducted through the condensate film gives the mass rate 
of condensation as: 


(rh) 


m 


A k-j 't ~ i 


n 


( 2 . 10 ) 


T W 
c> 


As the run-off is assumed to be gravity-dominated, viscous 
and surface tension effects being negligible, Bernoulli’s 
equation may be applied. It can be shown that the mass flow 
rate through a differential area P-dy is proportional to the 
square root of the local height above the plate surface y, . 

P being the perimeter of the plate allowing the flow of' the 
condensate (refer Pig. 1). Therefore, 


(dAl) out 


Ps|2gy . P d;v 


( 2 . 11 ) 


The total flow rate is found by integration to be: 

w) , 

47 • oy 
aw)/ 2 



(m) 

v ' our 

y = o 
= p iza -'f 

Jj = 0 

or. 

( m ) 

v OUT 

= -f- JIT ps 


( 2 . 12 ) 


This is similar to the inviscid result for the flow over a 


wier and 


has been obtained by Hollendorf and Chu [24]. 



(2.13) 


Considering the mass conservation, we write, 


= PA 


d 6 ( i ) 


QT 


(m) . - (m) , 

v 'm oub 

which, after use of has. (2.10) and (2.12), gives. 


, k l_ 
^ g 


(t q - t ± ) 

--cTtT*' 


% vi 2g ~f- [ & ( T ) ] 

'J • { 1 


5/? = a|ii). (? . 14) 


C'T 


Hence, using the dimensionless terms, from hqs. (2.y) through 
(2.8), the above differential equation may be written as; 

( 2 . 15 ) 


XX 


v X5/ 2 


0 


with the initial condition 

X = o at £ = 0 


( 2 . 16 ) 


where 


X = 


a x 

eft 


and, Y is a dimenoionle' s parameter for flow under gravity 

2 2 V 0 t/4 

•o .Ti- -•I u s ‘ 


5 ,nd is given by 


Y 


5 ^ A a v a P»-i fo , 


(2.17) 


It must bo noted that . the Y onl ^ T3aramet - 

that arises for f. low under gravity , and y = 0 represents bhe 
physical situation of no run-off. As y is a Strong function 
of P/A, whose value can vary greatly, y can take a wide 
range of values. It will be equal to zero, for a physical 
circumstance, only when P is zero. It can, however, ncmev- 



very small values, depending on fluid properties, the values 


of b, 0 g , cc, etc. This is an important point which has not 
been brought out by uollendorf and Chu [24] in their discussion 
on the flow parameter 'c:', which is very similar to y in the 
present case. 


For the no flow condition, i.e., y = 
equation (2.15), with the initial condition, 
results in the same equation as that from Sq 


0 , t he differential 
Fq. (2.15), 

. (2.5) for r = 0. 


It should bo noted that both r = 0 and y 


0 represent the 


no flow condition. 


As the expressions for heat flux for the transient 
f lira condensation on an isothermal surface are similar to 

those for other cases, to be discussed • later, they hove beer 
presented separately In section (2.1.6). 


2 . 1.2 


Plate^ with. Ij^sulat ec1_ 3o ttom Su rface^ and -. a Negligi ble 
Te):ap]e ratgire ^^a^oh’ XcT^ss] its. Thickness i 


The physical situation being considered is shown 
schematically in Fig. 2. A horizontal plate of uniform 
thickness b with an insulated lower surface, is initially 
at a temperature t^. At time r = 0, the top surface is 
subjected to a saturated vapour which is at a constant tempe- 
rature t_. It is further assumed that the Biot Humber, 
s 

B- = hb /k is small enough no that the temperature gradients 
1 P 

in the plate are negligible and that it may be taken to be 
at a uniform temperature. Evidently, low values of the heat 
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transfer coefficient, h, and high values of the plate thermal 
conductivity, k , with sufficient ly small plate thickness, b, 
will result in low values of Si. This situation .-.occurs very 
often in practice and physically corresponds to a negligible 
internal resistance to heat flow. In plates, the error 


introduced, by the assumption that the temperature at any 
instant is uniform will be less than 5 percent when the 
internal resistance is less than 10 percent of uhe external 
surface resistance, i.e., when hi < 0.1 [29 j. I'c is also 



7 out 
' ' 


Fig. 2: Co-ordinate system for transient film 

condensation on a horizontal plate, _ with 
insulated bottom surface and a, negligible 
temperature variation across its thickness. 

assumed that the length and* the width • of cne plate aiS 
sufficiently larger than the thickness so that the temperature 
gradients near the edge of the plate may be no leci-ed. 

Now we will consider the above situation for tlie 
earlier mentioned two circumstances of run-off. First, 
considering a specified fraction r, of the condensate, 
flowing out, we get the mass rate of condensation as; 


(m) ln 


= PA 


d 5( g) 


dr 


/ (1-r) 


(2.18) 
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Now, equating the heat released due to condensation to the 
energy change in the plate, we obtain, 


P I1.0 


d <5 ( T ) -up 

U-r) “aV“ = P p b C p, 




Mil 


d T 


(2.19) 


Since the heat flux conducted through the condensate film is 
equal to the heat released in the condensation process, we 


get. 


P h 


"f\or 


~TT-r T 


do_(jl 


t(r) 


dr 


k l f 1 ] 


(2.20) 


Using excess temperature and dimensionless film bnickneso 
and time, defined earlier, in the above equation, wo obtain, 


d 


x (?) = atf U) 


( 2 . 21 ) 


where 6 , the dimensionless temperature of the plate, is given 
by: 

^ e X llJi (2.22) 


°s 


' c s“ ' c i 


and, p the dimensionless condensation parameter is, 


P 


P h f a Ir 


i 


(2.23) 


Here, p is the most significant parameter, which controls 
condensation in the present case and in all other cases, to 

be discussed later. 
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Integrating the differential equation (2.21) and 
using the initial conditions that 


> inc o u 


X = 0 and $ - 0 at £ = 0 

t i at r = 0, us get. 


(2.24) 


X = 


LV. 1 ’.). d 
$ 3 


(2.25) 


iinilarly, Eqs. (2.2C) and (2.2?) yield 


X = 


(l~r) 


(2.26) 


Integrating :;;q. (2.26) with the initial conditions, Eg. (2.24) 
and making use of Ec. ( 2 . 25 ), we get. 


ti + Xi 

(1 - 0 ) 0 ' ” = 1 


( 2 . 2 ?) 


How , we will consider the flow dominated hy gravity. 
The heat released by condensation is equated to the energy 
conducted through the condensate film. This gives the mass 
rate of condensation as : 


(m) 


A 1 




{x) 


in 


fg 


U 




] 


The introduction of Bag. (2.28) and (2.12) into 
results in: 

5/2 


ti 


1 - X X - Y X 


(2.28) 

(2.15) 

(2.25) 



A balance between the heat conducted through the condensate 


film to the energy change in the plate , yields 

dt(-c) 

1 r; 

o(JT 


“ t ( T ) 

k 3 [ -- c v 1 


*u r* 

= P n 0 


* p 


Cl T 


(2.30) 


In terras of dimensionless ■variables, Eqs. (2.29) and (2.30) 


gives 


XX-;- (X % y x3/2 + p) \ + p y $/ 2 = 0 (2.31) 


with the initial conditions, 


X = 0 at K = 0 (2.32) 

The value of X(0) is infinite, as can be seen from Eq.(2.26) 

lor the no flow condition, y = 0, Eqs. (2.31) and 

(2.29) may be integrated. The result is the same as that 
obtained from Eqs. (2.25) and (2.27) for r = 0. The relations 
for no flow are: 


and , 


x= T 

(1 - 0 ) 


0 

0 + fh 

a = i 


(2.33) 

(2.34) 


It should be noted that all the above results can be 
used for a plate with varying thickness along its length, 
provided that the top surface is horizontal and that at any 
place Bi is not greater than 0.1. The only modification 
required is that b must be taken as Y/A, where Y is the 


volume of the plate and A the area of the top surface. 



2.1.3 Soni-inlini;iGL ^91^9'iL 

limes' aiicT 'Thiel:., Plates -. 

In this section, the film condensation ove j. a 
horizontal plate of uniform thickness, initially at tempe- 
rature t ± and behaving as a semi-infinite tody, will be 
considered. It is assumed that the heat transfer exxecus 
have not penetrated beyond the penetration dep-ch A, 
tho plato to-iperaturc at and beyond A is t j _, the initial 
t er.ro o rat uro. The plato is of a finite thickness b ano. 


, , .... r , he'-'sveo as a s eni—inf inix e oooy. 

for sliorr ^mes, x « jv-^v c ... -* / • > . 


( : ' a > in I 



( 

r-^t" 


2 


■ • - -ohem-tic of the transient film condensate n 

P,.,. times, on a horizontal plate* 0->- 

unixorm thickness and behaving as a s,mi- 
infinite body. 

Th o oni-infinite heat conduction problem i» solved by the 

. „ Aflr . .» ni „ no-ai rear heat transfer developed 
use of integral mo mods x 01 nouimocu 

_ T T97 2°1 A third order polynomial is employed 

by Goodman 1.27, 2oj. u 

for tho temperature distribution in the plate. 

Per determining the four constants of a cubic profile , 


given as : 


+ a „ s + a A 


3 


the following 


= a x + a 2 ~ + 3 

cond it ions are employed . 
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$> (a U), 5) = o 


(2.35a) 


agf 

ax 


(A , 5 ) = 0 


| 0 ; ( C. 


4(1, 5) 

as 


f (gL, € ) 


0 


( 2 . 3513 ) 


(2.55c) 


(2.55c 3 ) 


where X and A a: 


dime no ionless distance from, the top 


surface and the dimensionless penetration depth, respectively. 
They have been non-d imons ionalized villi, b as the characteristic 
C. ime ns io n . Xh e r of o re , 


\r 

JX 


mi! ■ 


A 


x/b 

A/b 


(2.36) 

(2.37) 


gf ( g) is the dimensionless temperature of the upper surface 
of the plate and f (# , £ ) is the heat flux, divided by 
thermal conductivity of the plate material. 

The first three conditions (2.35a,b,c) , known as 
natural conditions, arise from the considerations of the 
penetration of thermal effects in the plate and the boundary 
conditions at the top surface. The fourth condition is 
known as f-moothenina. condition and is obtained in the following 
manner. Conversion of heat conduction . equation s 


3 1 
3 t 


*2,. 

= a ■ -4 

3 x 


( 2 . 38 ) 



in terms of non-climensional variables, vields 


3 $ 


3^ 

air 


(2.39) 


•Ton, by differentiating Eg.. (2.35a) with respect to 5 and 
applying Ec. (2.39 ) * the condition given in Eg.. (2.35d) 
can be arrived [23], 

The appropriate cubic temperature profile is, 
therefore, obtained as: 


b = 6 0 [ i - ~ ] 

A 


(2.40) 


whore 


A = 


3 0 O ( 0 

z{i 0 & ) 


(2.41) 


Integrating heat conduction equation (2.32) with respect to 
X, as clone by Goodman [27, 28], we obtain : 


or 


f (0 O , 5) 


.4. f 2 
"”3 J ‘ 




[ 


'J 

\ l 


4, 2 w) , 

(t 0 . E) 


2 $ 


K. _M 

o d 5 f d 5 


(2.42: 


(2.43) 


How, first, wo will consider a specified fraction 
of the condensate flowing out of the edges of the plate. 

A balance between. the heat flux being conducted through the 
condensate film and the heat released in condensation of 
results in" 


the vapour. 
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P- a 


do(x). 


(l-r) ~dT 


V 


“l 


t., - t (t) 

r & O v ‘ I 

5 ( x ) ■* 


(2.44) 


1'Ion-dimensionalization yields an equation of the form 2 


' o 


= 1 


XX 


(2.45) 


As the heat flux conducted into the plate is same as the 
heat released in condensation, we obtain, 


"hjL. asjn. 

(l-r) cl x 


K H- (0, x) 

±' O 


(2.46) 


which, in its non-dimensional form, gives 


ty; (0 ’ 


'(I- rj X 


(2.47) 


'from Ihin. (2.55c) and (2.47), wo get 


(2.40) 


The introduction of Sqs. (2.45) and (2.48) into the T5q.(2.43) 

yie3-dn ,1 


{ (l-r)'” - ^ X"] X + C 2 (l-r) 


4 


1 


2 XX j X J = 0 

(2.49) 


P r0 ,-, -bi'ie above equation it can be shes®. that XX must 

* 

be a constant, since r and p are constants and X and X 
vary with tine. 


XX = Constant = % 


2 . 50 ) 
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then from 33a. (2.49) 


X = 


Q X- 


where 


o = : r ] 


(2.51) 


(2.52) 


(1-r)^ - %“ 

Integrating the differential equation (2.51) and using the 
initial condition, Eq. (2.52), which io valid in the presonl 


case also, we get, 

>• JTR 


Hence 


TJ 

i! i 


l 


(2.55) 


(2.54) 


Using Eq s . (2.52) and (2.54), we get the quadratic equation 


2 - [2(1-1-) + 4 P 2 ] % + 0 - - r ) 2 


0 


(2.55) 


% - L ' ' 1 ^ 

whose roots are given by 


T\L 


(1-r) 


6 2 + 




4 0/1 /I. 

I (l-r) + Of P 


( 2.56 ) 


Hence, the solution for Sq. (2.50) and, thus, xor 
the differential equation (2.49) is given by, 


X = H K 


1/2 


(2.57) 


where 


T'T 


= Jiih = [2 (l-r) + 3 P ' + 


'1 ~ l- - V - . - ' ^ 1 — N 


16 (l-r) p 2 + 4 P 4 ] 


1 / 




9 

( 2 . 58 ) 
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Since, the rate of increase in the condensate film thickness 
should not he more than that over an isothermal surface at 
the initial temperature t ^ , because of the plate thermal 
resistance introduced in the present case, re obtain tne 
following condition from Sq . (2.9), 


H < |2(l-r) 

Hence, the accepted value for rl is. 


IT 


r , N 4 2 r 3_S f i N c.2 ~ IlK r4 

r 2(l--r) + ^ P ~ \i r K~ ( x ) >' 9 ^ 

J (2.59) 


Prom the Eqs. (2.41), (2.45), (2.48) and (2.57) the non 
dimensional penetration depth is obtained ass 


A = M 5 


:l/ 


(2.60) 


where , 


xi — pT-T 


Eqs. (2. ;5) end (2.57) -give the dimensionless 
the top surface of the plate as ; 


<b. 


IT 


2 


2 ( 1-r) 


(2.61) 
erature of 

(2.62) 


Prom Eqo. 


(2.40), (2.60) and (2.62) the temperature in the - 

plate at. any distance X is obtained ass 


$ = r 1 - — ] (1 

2(l~r) 


n c -i/'2)3 
n ' 


(2.63) 
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Ths equation (2.62) is quite significant and shows 
\1m\# w 'O *■ *■ l_) ore, "cure of tho top surface ox uho pla^e will 
attain r constant val_ue immediately after it is exposed to 
the saturated vapour , till the plate behaves as a semi- 
infinite booy. This tennoe nature is a j. unction o.l r cmd p. 
For the no flo : condition, it is a function ox p only. 

'.‘or the flow dominated by gravity, the mass rate of 
oond onset ion is obtained ass 


(m). 


' u o v M i 

nr-cT J 


(2.64) 


Iron ten. (2.12) .and (2.13), which are applicable in the 
present case also, and Bo. (2.64) 0 O iB obtaino<3 111 terms 


of X and y as; 


1 -- XX 


(2.65) 


Again, equating the heat flux concocted through me conden- 
sate film to the heat conducted ih" c ° the pis. uo at ^ - 0, 


gox , 


t 0 ^> 


. 2J. (0, t) 


( 2 . 66 ) 


... -f-PTins of dimensionless 

hro res sing the apove equation m c® ruo 


variables and employing Fas. (2.65) ano ( 2 ‘->oc) , -re t 


f($ c , % ) = p ( X + y X 5/2 ) 


(2.67) 


The introduction of Ben. (2.65) ^ n ° into ' ciie 

which is basical l; r a conduction equation xor. tne plat^ 
therefore , applicable here also, results 



31 


[1 - (XX + Y X 5/2 ) 2 ]X + | f (x + Y ^ /C ) 


4 „2 , .. ,3/2 ^ 


(1 - xx- Y x' /2 ) (2r + — Y 

+ l y 2 X 5 X + | y X l/2 X ) = 0 


11 .. x 3/2 ;2 


( 2 . 68 ) 


with the initial condition given by Ecu (2.32). The penetra 
tion depth is, thus, given by Egs. (2.41)? (2.65) and (2.o7) a 


T _ _3. r .... 

A ~~ ~ 


P L (x +7?^’ ’ 


X ] 


(2.69) 


The temperature at any distance X is obtainec from nos. (2.40) 
and (2.69) as : 


fl (X-l-Y X^ /2 ) 

$ = 0 O [ i - Hr 7 : 5 T 2 T x j 

0 ' (i - xx- y \ ) 


5 


(2.70! 


Again, for the no flow condition, y = 0, Ecu (2.68) 
yields the same expression as that obtained from Ecu (2.4./) 
for r = 0. 

2.1.4 Plate of Fi ui te .ThAofeiggg an d Bo ttom _Surf .aco^jit 
Pi'xc d* Temp erat ure ; 

The transient film • condensation on a horizontal 
plate? with finite but uniform thickness , and whose surface 
in maintained at the initial temperature t ± is considered 
in this section. The solution of thin problem will follow 
the previous solution, in which the plate had been considered 
as semi-infinite. Hence, the finite thickness of the plate 
will affect the solution beyond a time interval given by 
the time at which the penetration depth equals the depth b 
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for the semi-infinite solution. The physical situation is 



Again, for the present problem, the integral method 
developed by Goodman [27], is employed. A third order poly- 
nomial for the temperature distribution in the plate is 
assumed.- For determining this cubic profile, of the tempera- 
ture in the plate, the required four conditions are; 


0 (0, C) = 0 O ( £ ) 


(2.71a) 


| 4 (0,5 ) = - f (0 O , 5) 

O 4s. y 

<t ( 1,5 ) = 0 
M ( 1,5 )= 0 

«\ **r 


(2.71b) 

(2.71c) 

(2.71d) 


The first three 


conditions are derived from the boundary 


conditions. As the temperature of the 
fixed and does not change vith respect 


bottom surface is 
to tine, we obtain. 


kA 


2=1 


0 


(2.72) 



and, hence, from the heat conduction equation (2,39), the 
fourth condition, given by Hq. (2.71c! ), is obtained. 

The appropriate cubic temperature profile is, therefore, 
given by [27], as: 

0CK, ? ) = [ | 0 O (?) - | f(0 o , £)] (1 - X) 

+ \ [f(0 o , o - <6 0 ] (1 ■■ X) 3 (2.73) 

Integrating the heat conduction equation (2.39) with respect 


to X and applying the necessary conditions as done by 


Croodman [27], we obtain, 

12 [f (0 o , K) - 0f o (5)] 


cl0 Q (?) 


df (0 O , C ) 

“ “ n 


(2.74) 


Tor the first type of flow, i.e., a specified fraction 
of the condensate flowing out across the edges of the plate, 
the Bqs. (2.45) and (2.48) are applicable, and an introduction 
of these equations into Eq. (2.74) results ins 

(p -l- 5 X) X -i- (5 X + 12 X + 12 p) X - I2(l-r) = 0 (2.75) 


The initial conditions for the above equation are to be 
obtained, from the semi-infinite solution, Eq. (2.57). Values of 
1 and X at the end of the semi-infinite solution, i.e., for 
x = % ,, when A - b, are used as the initial values for the 

above equation. Hence the ' initial cond it ions are j 



at C = ?o 

w 


X 


when A = 1 (2.76) 
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The temperature at any distance X is given by Eq L s. (2.45), 
(2.43) and (2.73) as; 

C { ( 1-r ) - XX} (1 x - 2~) 

-,2 

~ P X (X - ™)] (2.77) 

For flow under gravity, Eqs. (2, 65 ), and (2.67), of 
section (2.1.3), are applicable . Introducing these equations 
into Eq. (2.74), we obtain, 

(p + 5 X ) X + (5 X + %• y \^ 2 12 X-i- | p y X ly/2 

12 p) X - 1 - 12 y (P +X ) X 5 ^ 2 ~ 12 = 0 

(2.73) 


with the initial conditions given by Eq. (2.76). The values 
of , X„ at in the present case, are obtained from 

Eq. (2.68). For the no flow condition, Eqs. (2.78) and ( 2 . 75 ) 
are identical. The temperature in the plate is obtained by 
the substitution of Eqs. (2.65) and (2.67) into Eq. (2.73) as; 

0 = (1 - : .) [(1 ** XX - Y X 5/2 ) (1 + X ~ “2') 

- p ( X+ Y * 3/2 ) (X - ^ )] (2.79) 


2 . 1.5 Plate of F init e__ Till c igies_s_ and 

The physical situation to bo considered is as shown 
in Fig. 5 . A horizontal plate of uniform thickness b, with 
an insulated lower surface is initially at a temperature t ± . 
At time x = 0, the top surface is exposed to a saturated 
vapour, maintained at a constant temperature t g . Initially, 
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the plate behaves as a seni--ia.tiaite body in the x-direction. 
After a time n, when tne ponetr. ... tion depfh A ecrud.. ,o , 


finite thickness of the plauo will affect tne oducion. A 
solution will bo developed for time x >x g . This problem 
would also arise when a pi?/* of finite thickness is immersed 
in a vapour so that condensation occurs on both sides, the 

adiabatic for symmetric heat transfer conditions 


mid plane being 



(m) 


out 


. p. r i .« - n «<!•!• cm and physical model xor che 

Pi,v. 5: Go--oraiiiatc d - 1AU . 0 - 0 ipte of 

transient file condensation on c, pnave o 

Ukj X ^ -I ** ^ ~ . . n ^ 4.0 ^ oo Q coni BU.L Xd.cc * 

finite thickness and 

. -i -i -rin'r^ i integral mo tlio o ? 

Again, for the present problem, uhe mte 0 i 

Ti.-.ooT' he^t transfer problem 

discussed earlier, for the non-linec.r he.. 


JL 
„ 0 

UIGUULPW'-’U -a- y — 

is employed. As only three boundary conditions for me 
■tevipG3.-n.turo in tie plate are inorn, namely 

0 (0, 5) = 2S 0 (5) 

U <0, E) = -f(0 o »s) 

a ^ 


( 2 . 80 a) 

( 2 . 80 b) 


\& (1, 5) = 0 


(2.80c) 

¥e can Imploy only a second order profile- The appropriate 
quadratic temperature profile for che plate 
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,2 


0 = [ 0 O (S) - ?r f (0 Q f C )] + | f (0 Q , 0(1- I)' 

(2.31) 

Integrating the heat conduction equation (2.39) with respect 
to ji. and appl'ging the required conditions, we get, 


ri cj 


80 
3 


If (1, 5) -|f (0,5 ) 


which, after substitution of Eos* (2,80c) and (2 # 8l), results 
ins 

d r f ^ \ 1 . \ i (282) 


or 


£0*o’ e } = cfT C^oU) - I 5 )] 


*«L. X «<*<..<> 


(2.83) 


o' * """cTs “3 clT 
The Eq. (2.02) is basically a heat conduction equa- 
tion and can 'be used for the purposes other than the present 
ones, for example, for nonlinear heat transfer, in a slab of 
finite thickness and one side insulated, with the temperature 
of the -other surface a function of time and the heat flux a 
prescribed function of surface temperature and time. 

■■or the flow, situation , when a specified fraction of 
the condensate is flowing out across the edges of the plume, 
Eqs. (2.45) and (2.48) are applicable and use of these equa- 
tions alongwith Eq. (2.83) results Ins 


(X + —M- ) X + ( X+ j3 ) 


X = 


0 


(2.84) 
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with the initial conditions given by Eq. (2.7 6). The values 
of X s and X f , at £ = K a are to be obtained from Eq.(2.57). 

Apparently, thq above equation does not involve the. 
flow terra r. But the value of r affects the solution of this 
equation, because the solution has to follow the semi- 
infinite results, which depend on r, as shown earlier. The 
dimensionless temperature is given by Eq. (2.81) after 
employing Eqs. (2.45) and (2,48) as; 

v-2 : 

0 = 1 - [ X - p (I - -y )] ----- (2.35) 

(1-r) 

Por the flow under gravity, Sqs. (2.65) ancl (2.67) 
of section (2.1.3), arc, again, applicable . Hence these 
equations with Eq. (2.03) gives, 

(# + X ) X + (X + | Y X^ 2 + ~ p Y x 1 ^ 2 + j3) X 

+ p y x 5/2 = 0 (2.86) 

with the initial conditions given by Eq. (7.16). Here, the 
values of X and X 0 at £ = are given by Eq. (2.68). 

Prom Eqs. (2.65)* (2.67) and (2.81), wo get the dimensionless 
temperature in the plat s as ; 

= 1 - [ X + f, (2 - K ) ] ( X+ Y X 5/2 ) (2.87) 


0 
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2.1.6 Heat Flux: 


Since the equations of heat transfer, for the transient 
film condensation on a horizontal "body, with various boun- 
dary conditions discussed earlier, are similar, except for a 
change due to the type of flow, they can be presented at one 
place. 

For a specified fraction of the rate of condensation, 
flowing out across the edges of the condensing surface, the 
heat flux released due to condensation of the vapour is given 
bv: 


P h 


fg 


*!■ v ■ v 

Vj'-rf "• 


( 2 . 88 ) 


which, in non-dimensional form, gives. 


v a 

^ rx 

q = - ' F 1 • 


X l~rj * x 


(2. as) 


ho may, for simplification, introduce a heat flux parameter, 
q given as : 


Q. = 


<V B - ) 


(1-rl A 


( 2 . 90 ) 


The above equations are va,lid for all the five different 
boundary conditions, discussed earlier. 

For flow unde:; gravity, the heat xlux conducted 


through the condensate film is given by, 

r 1 

q = % [ -- g-:- j 


( 2 . 91 ) 
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which, in terms o f ciim aim ionless variables, gives . 
q = -j 3- P (i + Y X 3/2 ) 


( 2 . 32 ) 


Hence, the heat flux parameter q is obtained as: 

1 = = ,3 ( A-:- Y A 3/2 ) (2.93) 

r “9 s > 

( ~- v ~) 

It should be noted that the top surface temperature t Q 
in liq. (2.91) is equal to t^ in case of condensation on 
isothermal surface and to t(r) for condensation over a 
plate with bottom surface insulated and temperature variation 
negligible . Again, in this case, the above equations are 
valid for the various bound ary conditions, for the plate, 
discussed earlier, provided the flow is gravity dominated. The 
heat flux parameter’ q, for the no flow condition, is, therefore, 
given by, 


A 


(2.94) 






Dimensionless plate temperature 





Dependence of X on y for transient film - condensation 
on a horizontal plate, with bottom surface insulated and 
temperature variation across its thickness negligible. 




2 • 2 Be s uit s_ and Discus s ion ; 


(a) 


of 


Isothermal Surface : figure S 
the dimensionless film thickness 


presents the variation 
X with respect to 


dimensionless time £ , for various values of r and y. Haile 
for any value of r, which must obvious ly be less than 1.0, 
the film thickness X increases with time, it achieves a 
constant value after a certain time for y >’ 0. The time, for 
X to become constant, is found to decrease as the value of 
y increases. For the no flow condition, r = 0, Fq. (2.3) is 
the same as the expression obtained bj* Gebhart [19], and by 
P.olis enow and Hartnett [23b]. For small values of C A t/h.* 

J? 1 -1- o 

the results given by 2q. (2.3), or by Eq. (2.9), are in good 

agreement with the exact results obtained by Pfahl [22] for 

the no flow condition. The error is only 1.5 percent when 

CL A t/h~ eouals 0.135. The equations (2.5) and (2.9) are 
P ]_ -*-6 

also much simpler than the expressions given by Pfahl [22], 

(b) Plate with Insulated Bottom Surface and Negligible 
Temperature Variation Across Its Thickness; 

The calculated dimensionless film thickness X and the 
dimensionless plate temperature 4> are shown in Figs. 7-9, 
for a range of r and y. Increasing r or y is seen to decrease 
the maximum film thickness and to increase the rate of tempe- 
rature rise of the plate. This is reasonable since a thinner 
film is expected for greater flow, obviously, offering less 
resistance to heat flow. 


In each ca.se the plate temperature 




■ ' 


1 — 

- JU' 




IIS 


a horizon 
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reaches a value very close to saturation temperature of the 
vapour. While for any value of r, the film thickness is 
constant after a certain time, it becomes zero ultimately 
for y > 0. This is expected since no condensation occurs 
when the body reaches the saturation temperature. Clearly, 
therefore, physically, the model based on gravity is more 
suitable for large times, since the model based on flow 
fraction r gives a constant value of A . At low times, this 
latter model may be quite suitable, particularly, in practical 
situations. Mollendorf and Chu [24] have also obtained the 
expressions for the film thickness and the plate temperature 
and have presented the results for various values of the flaw- 
parameter 'of, which is similar to y in our case. Their results 
are same as ours. They have also reported that for the no flow- 
condition, the agreement between theory and experiment is 
excellent . 

Figure 10 shows the effect of j3, the dimensionless 
condensation parameter, on the film thickness for the no flow 
condition. The rate of heat transfer increases with increase 
in p, and, hence, the plate temperature t reaches a value 
close to the saturation temperature t g much earlier. Obviously, 
the maximum film thickness is greatly reduced as (3 is increased. 
The same trend is obtained for all other cases, i.e., various 
boundary conditions assumed for the plate, and for all values 
of r and y. 




Variation of X with f, over a range c 
for transient film-condensation on a 
horizontal plate of uniform thickness 
behaving as a semi- infinite body, for 
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(c) Solution for Short Times and Thick Plates : 

Pith. an assumed third order temperature profile in 
the plate, the semi -infinite solution, for a fraction r of the 
condensate flowing out across the edges of the plate, Bq.(2.57) 
indicates that the film thickness is proportional to the 
square root of time, the proportionality constant being a. 
function of r and f., i.e. , of the fluid and the plate proper- 
ties, At, flow fraction r etc. The constant of proportionality 
IT, given by Eq. (2.59) decreases with an increase in r and p. 
Pig. 13. shows the film thickness for a range of r. Bor the no 
flow condition, Pfahl [22] has reported the same behaviour 
for the film thickness. There is a difference in the propor- 
tionalitjr constant, the evaluation of which is quite difficult 
in his case and involves the error function. For small values 
of CL At/h. P „ the results given by Eq. (2.57) are in good 

- L O 

agreement with the exact results obtained by Pfahl [22]. For 
a saturated steam at atmospheric pressure condensing on a 
Hi-steel plate (Hi - 40?:'), initially at 27°0 , the results 
obtained in this study differs by only 0,1 percent from those 
obtained by Pfahl [22], 

Interestingly, the temperature of the upper surface 
of the plate $ , given by Eq. (2.62), attains a constant- 
value, depending upon r and j3, as soon as it is exposed to a 
saturated vapour; and remains at that temperature as long as 
the effect of the bottom surface, on the solution, is not felt. 








mensionless penetration dealt 



For all value 
of r and j3 


dence of d 1 rn e ns io nS ess penel r at i o n 
r a nd y f nr t if n n s sent f I ! r?i ~ c on de 
horizont ! , ^Having as a se 



The greater the values of r and j3, the larger is the surface 
temperature 0 Q . 

For the flow dominated "by gravity, the effect of y 
over the film thickness is shorn in Pig. 12. Here, the surface 
temperature 0 is not a constant but increases with time. 
Greater the value of p or y, greater is the rate of increase 
of 0 O . 

Por the no flow condition and a fraction r of the con- 
densate flowing out across the edges of the plate, the pene- 
tration depth, given by Eq. (2.60), is proportional to the 
square root of time, as may be expected. The constant M in 
Eq. (2.60) is a function of r and j3 and is found to be the 
same for almost all combinations of r and j3 , and, hence, the 
penetration depth, for all values of r and j3, 'is represented 
by only one curve, shown in Pig. Ip, As the temperature of 
the upper surface of the plate t , in this case, attains a 
constant value as soon as it is exposed to the condensing 
vapour, it is a case of a step change from the initial tempera- 
ture t^ to a constant temperature t 0 . Using a third order 
temperature profile, Eckort and Drake [31] give the dimensionle; 
penetration depth A as >jl2 \ for the above case. Goodman [27] 
has used a third order profile of a different form and has 
found a as equal to 424 k * which is expected to bo a better 
result than the result obtained by Eckert and Drake due to the 
assumptions made. In the present study, M is within 1 percent 

of J 24 for. almost all combinations of r and )3. Dig. 13 also 

N j,!. 




'■ ' 
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shows the effect of the flow parameter y on the penetration 
depth. It is observed that the rate of penetration decreases 
as y increases. 

(d) Plate of finite Thickness and Bottom Surface at a. Fixed 
Temperature ; 

Fig. 14 shows the effect of r on the film thickness , 
which- is found to continue increasing with time. It is also 
observed that the switch over, from semi-infinite solution to 
the solution for finite thickness , is quite smooth. The devia- 
tion from the former is a gradual process. Fig. 15 presents the 
temperatures at different positions across the plate for the 
no flow condition. The upper surface temperature, which is 
constant as long as the plate behaves as a semi-infinite body, 
starts decreasing as soon as heat transfer occurs at the bottom 
surface , and after a very short time, starts increasing again 
till a. maximum value is reached. Temperatures, at other 
locations considered, also show an increase with time till 
the maximum values are attained, and , then, indicate a 
decrease as for the surface temperature. The decrease in tem- 
perature is expected from the decrease in heat transfer to 
the body as the film thickness increases. The energy lost 
at the bottom becomes larger than that gained at the upper 


surface . Obviously, t(x) is greater tha: 


■a. i or x < o. 


5ven 


after a very large time. The exactly similar behaviour is 
observed for any other value of r. 




• '■■■ ' »■ .‘'J| 
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SSEBf®; 




X for transient film-condensation 
finite thickness and bottom surface 
berature 


Effect 
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For the gravity dominated flow, the film thickness 
attains a. maximum value after a certain time and then 
decreases to a constant value which is maintained for all 
times. This steady state indicates that whatever vapour is 
condensed, flows out and the heat gained at the upper surface 
is lost at the bottom surface. As y increases the maximum 
film thickness decreases, as shown in Fig. 16. Similarly, 
the temperatures at different locations in the plate, for 
y > 0, are constant after a certain time, as can be seen 
from Fig. 15. A comparison has also been ma.de between the 
film thickness solutions for a semi -infinite plate and those 
for a plate of finite thickness, as presented in Fig. 17. 
Clearly, the difference is small over the time period 
considered. 

(e) Plate of Finite Thickness and Bottom Surface Insulated: 

The solutions for this problem are an improvement over 
the solutions for the plate with insulated bottom surface and 
negligible temperature variation across its thickness. The 
behaviour of the physical variables is the same except for a 
small difference in the values, resulting from the assumed 
quad rati c temperature profile in the plate. Fig. 13 shows the 
effect of r on the surface temperature and Figure 19 
presents the behaviour of the film thickness for a range of 
y. As can be seen, the temperature approaches the saturation 
vapour temperature and the film thickness gradually goes to 
zero for the gravity dominated flows. 















= q / ( kp-0 s / b) 
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A comparison has been made between the results obtained 
for the condensation on a plate, of uniform thickness and 
insulated! bottom surface, with and without a temperature 
variation across its thickness, as presented in Pig. 20. The 
behaviour of the film thickness is quite similar in the two 
cases, though the temperature variation in the plate is 
obviously quite different. 

(f) Heat Flux: 

Fig. 21 shows the variation in heat transfer rates, 
with time, for various boundary conditions of the plate, for 
the no flow condition. Heat transfer rates are found to 
decrease sharply with time for all the cases. For the no 
flow condition, the film thickness increases with time and 
this would obviously result in the reduction of heat transfer 
ra.tes. This result would be applicable for short times when 
flow has not started due to surface tension effects and for 
circumstances where a wall obstructs the condensate flow. This 
figure simply shows the similarity in the heat transfer 
results for the various cases. Similar curves may be obtained 
for other circumstances. 


XAIIIiIA.il PIIH C 0 ITDEUS AT I Oil FLOP ON A HORIZONTAL 
ISOTHERMAL SURFACE 


3.1 Analysis,; 

The present chapter discusses the analysis and the 
numerical results for laminar film condensation, flow on a 
horizontal isothermal surface (Fig. 22), subject to the 
approximations discussed in section (2.1). The analysis 
ignores any disturbances, caused by the condensation of 





o 




min 


?ig # 22 Physical model and co-ordinate system for 
laminar film condensation flow. 

vapour, on the velocity and temperature . of the condensate. 

The complete analysis has been divided into three parts, as 

already discussed in Chapter I. 

3,1.1 Transient Pila. Co nd ®I*sat io P -. ; 

When a horizontal isothermal surface, at a temperature 
t-, is exposed to a saturated vapour, at a uniform temperatur 
t corresponding to its pressure, condensation occurs. The 

s 
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liquid , thus condensed, does not flow out across the edges 
of the condensing surface unless the condensate film is 
sufficiently thick, so the t the hydrostatic forces overcome 
the surface tension and viscous forces. Hence, transient 
film condensation, with no flow occurs initially and is 
considered first , as given below. 


The thickness of the condensate film, 5, which is 
aero at time x = 0, will continue to increase till the flow 
starts. Therefore, from Eq. (2.3), the condensate film 
thickness for no flow, is obtained as; 


\| 


2k l (^s " 

p h,;; T 


(3.1) 


X 


o 


If the condensing surface is terminated by a sharp 
edge, the film is essentially a large, two-dimensional 
drop [21], whose height increases with time. The shape of 
this drop, just before the flow starts, is determined by a 
balance between the gravity and surface tension forces. 

Bankoff [ 25 ] has analyzed the problem of spreading of large 
drops and both his analytical results and experimental 
observations show that the surface tension controlled film 
is of practically constant thickness over most of its length, 
when 6 . the film thickness at x = 0, is small as compared 

to the length of the plate. Por water with a contact angle of 
%/ 2 , for example , the film thickness attains 99 percent oj: 
its maximum value within 3 <S maz of the edge [21]. 



0 . 1 - 
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The expression which determines the film profile 
23) has been obtained by Tankoff [ 25 ] and is given by. 


= 3 2 log \ 


( x |l - + 1)] 


Constant 


(3.2) 


whore , 


(3.3) 


( P- P V )S 


and ; 0 is the surface tension. The sign in front 01 uhe oq.ua ie 
root terms would be reversed for a drop advancing to the ^ 
right. 3? or the co-ordinate system with x = 0, when y = 

tho coa-tesrt in Bo.. (3.2) [25]- K5- 2? ? reselrts 

profile of a. largo drop of watey, at different toupoi.-tuies, 

advancing to the left. 

o) w iU give the 

If end effects are ignore a , nq. U> • 2 J wi-j- o 
ftln thickness a,, at * = 0, which must be attained before the 
flow occurs. AS the variation in the film thickness at the end 
of the condensing surface is very small, as discussso earlier, 
we may obtain the time t 1 required for the film thic-nes,, to 
bo Cb, assuming it to be uniform over the entire surface, 

from bq. (.'#!)* 

3.1.2 T runs i ent > laminai^ film .Cpndej^a_oipn_ |'lo}r •> 

After the flow of the condensate has just started, 

Ty, ihio section, the 

it would undergo a transienvproceau. In uai, 

„ TT wpen considered, 

transient, laminar f ilm condensation xlow hep . 
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Assuming the flow to be unidirectional, a simple balance 
between the pressure gradient and the viscous and surface 
tension forces, at any instant, gives, 



where, u is the 
due to surface 
sate fluid arid 




0 


(3.4) 


velocitjr in the x-clirection, p a the pressure 
tension forces, jj. the viscosity of the conden- 
5 a function of x and t. 


Since the thickness of the film is assumed to change 
substantially only in the direction of flow, i.e., in the 
x direction, p a is given by [ 32 , 33] ; 


P 


a 


a 


*7 


(3.5) 


where, is the radius of curvature of the film surface and 
is related to its thickness by the well-known equation [32 ] % 

r,, = - [i + (|4) 2 f /2 / ^ < 3 - 6) 


As << 1, we may assume that 


P 


a 


3 2 6 
cj — ’2’ 
9X 


(3.7) 


Hence Eq. (3.4) may be written as; 


.2 

n 9 U 

9 y 


f n — P ) v + o — v = 0 
IP V ° 3x 9 x 3 


- 


(3.8) 


As the velocity at y = 0 is zero (Fig. 22), and the interfacial 
shear between the liquid film and the vapour has been assumed 

to be negligible, we get, 
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u = 0 at y = 0 


(3.9a) 


3 u 
3~~y 


0 at y 


(3.9b) 


Integrating Eq. (3.8), over y, and using the above boundary 
conditions, we obtain. 


u 


.1 

n_ 


(p 


3 _5 
3 x 


D o 

g ~~ - a ) C J ~2 ~ y<5) 


( 3 . 10 ) 


3 x" 


where , 


P = ( P~ P v } 


(3.11) 


For an assumed linear temperature profile in the 
condensate film, we write, 


ai 

3y 


' G i 


T(1T,tT 


A t 


( 3 . 12 ) 


where , 


At = t c 


Integrating Bq. (3.12) from y to 5, we get, 
t„ - t = (1 - ) At 


(3.13 


Considering the energy balance of an element of the 
film of height 6 and thickness dx, depth being unity, we 
write, 

6 


IE w _ u 
3 x ‘ n fg 1 3 y 


y = o 


P C t._ / (t - t g ) dy 


Pi 


+ r| J u (t 
3x 0 


t ) dy] (3.14) 
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where , 


3 m 

af 


; Condensation rate per unit length. 

Using Eqs. (3.10), (3.12) and (3.13), the above equation may 


be written as 
3_m _ 


1 p 
^- [ 7r - P 0.. U + — P 


h 


P-j 9 T 


!i~ 3- 


{ 6^ (pg - or 


3 


h" jj lu 32 


a x 

(3.13) 

Now, considering the mass conservation for the element 
we obtain, 

6 

1 

0 

Introduction of Bq. (3.10) into the above equation yields 


3 — f)>] 


8m . 8 6 i 

h = [ l p u dy + PaV ] 


(3.16) 


3 a 35 IP 3 

a x ~ p f \i 3 2 

3t 




(3.17) 


a x' 


From Eqs . (3.15) and (3.17), we get 


p ( h fv + \ °v, At) 


fg 

Tc7*'A’t* 


9_6 p 
9 x ~ 3(i 

,3 > 3 


(h f „ + 4 c „ At) 


pi 


h At 


At 

1 P 1 


< 6 3 ( p g - cr “3 )> ] = 1 (3.18) 

u o 2C 


or, 


V . (1 + USs.J. 6 + 6 3 [ -- { 6 + 3 |4 } 

pg L pg 


a x 


{ 6 + 3 (||)^}] = 


. c o 3 vh, At 

9b\2 -i-i _ ± 


3x 


p g K 


•fg 


(3.19) 



which, in terms of dimensionless distance, X, yields. 


C At 

3P v (i + 1 ^ ) - 2 


3^6 


P g 


"i'g , 3 5 , ,3 r a { c 

° ST + 0 l 6 “4 

p gL o ii 


+ 3|l^l}~{6^-|+3 (f|) 2 ] 

9 ll ^ *x r J *\ ~T*~‘ 3 ^ 

d ii a a 


3 v At 2 


p g h. 


fg 


( 3 . 20 ) 


whore 


c/3 


( 3 . 21 ) 


Hence, we obtain. 


R 56 + 5^ (66 1V + 36- 6 m ) - (56” + 36’ )] = EL 


r , 2 ^ J _ rrT 2 
( 3 . 22 ) 


whore 


C At 
1 p l 

3 p v (i + 1 


- ) 


P g 


K = 


3 v k^_ At 


P g kfv 

•Lg _ p 

p/yjj 

Bond Humber, Bo = 


(3.23) . 

(3.24) 

(3.25) 


a 



If surface tension effects arc ignored, or for large 
values of Bond number, Eq. (3.22) yields, 

R 6 6 - 6^ (66" + 3 6 ’ 2 ) = KL 2 (3.27) 

How, proper conditions must be determined to solve the 
above differential equations. For a condensing surface ter- 
minated by sharp edges, the initial condition is obtained 
from the section (3.1.1), i.e. , 

6 (X, 0) = 6 X (3.28) 

For other end conditions, i.e., for a rounded fall at the end 
of the condensing surface, we must obtain the film profile, 
necessary for the flow to start, by some other methods. For 
certain fall conditions, it may bo expected that the film 
thickness at the time of flow to start is quite small. 

How, since the flow is symmetric on both sides of the 
condensing surface, we get, 

U. (0, t) = 0 (3-29) 

The other boundary condition, required for the solution of 
Eq. (3.26), must be known at the end of the condensing 
surface. For certain fall conditions, we may take r) min as 


cons tout, where. 
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*0 


(3.31) 


The other possibility is to take 


' ( 1 , t) 


Constant 


(3.32) 

as a boundary condition, because the rounded fall at -Kh~ 

rao ond 

has a definite slope. Using the initial condition, pa .(3 29 ) 
and either of Eqs. (3.31) and (3.32), as applicable 5 xr 
o b t a i n th e t rans i o nt solution from Sq . (3.27). 


're may 


If surface tension effects are also considered, 


¥0 


must solve the equation ( 3 . 22 ), which requires two m oro 
boundary conditions to be known. As the flow is syminei ric 
on both sides of the condensing surface, the fluid a -j- 


= 0 


wi 11 be at so 


3 3a 


ro velocity, and, hence, from Eq. (3.1 q) ? ^ 


- f (0, a 


0 


(3.33) 


Therefore, for a fall condition, where, both p min and 5 '(i sT ) 
are constant, Sq L . ( 3 . 22 ) may be solved using the initial 
condition and Eqs. (3.30) through (3.33). 

Since, the temperature profile in the condensate film 
is assumed to be linear, the local heat transfer co-efficient 


h„ is given by; 

.A. 


x 


k. 


6Xx , ' T) 


h(x, t) 


(3.34) 
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And, hence , the local Musselt number is 


HU, 


rrxTxy 


(3.35) 


3.1. 

► 3 

Stead#*. Lamina: 



for 

negligible 

values 

of - 

j O , -jO « \ 0 . 

for 

st 

eacy 

f laminar f 



a 3 

[65" + 3(5' 


-KL" 


s ati on . Flow : 

on effects , or for large 
c ond ens ate film p rof ile 
on flow as? 

(3.36) 


with the conditions 
o’ (0) = 0 


y, = Constant 

'min 


(3.37) 

(3.33) 


If the transient solution is obtained, 6 max will be a known 
quantity, and as such, 

*(* o) - 6 (5.39) 

may be used as a boundary condition, instead of Eq. (3.38). 

Even, if the film thickness at any X is known, Eq . ( 3 . 3 6 ) may 

bo solved with the help ox Eq. (o»o7). 

Nimmo and leppert [20,21] have also analysed the 

steady, laminar film condensation flow, with negligible 

surface tension effect, and, have obtained on equation similai 

to Eq. (3.36)., Using the condition that the film thickness at 

Y n d is known and that given by Eq. (3.37), ihey have 

A - u > °max’ 








presented the exact solution of Eq. (3.36) in terms of 

elliptic integral. 
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If surface tension effects are considered, Bq.(3.22) 

yields , 

& [g~ (66 1V + 3 o' 6"') - (66" + 36' 2 )] = KL 2 (3.40) 
The above equation raay bo solved employing the conditions, 

6'" (0) = 0 (3.41) 
and 6' (1) = Constant (3.42) 


in combinatio n with the other two boundary conditions dis- 
cussed above. It may be noted that Eq. (3.36) or Eq, (3.40) 
may be directly derived from the equations of motion, mass 
conservation and energy balance. .Again, in the present case, 


the local heat transfer co-efficient is 



(3.43) 


and, the local Husselt number is 


3tfu, r 


X 

oTx) 


(3.44) 


3 . 2 Results and. Discus sio n % 

(a) Steady, Laminar Film Condensation Flow? 

The condensate film profile for steady, laminar film 
condensation flow, with negligible surface tension effects, 
is given by Eq. (3.36). Fig. 24 presents the numerically 
obtained film profile for various values of '0 m ^ n . There is 
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very little difference in the film profiles for t) ±r < 0.4 
'The sane behaviour has also been reported by STirnmo and 


Lepport [21]. In this range of r u ^ s 


i m i n ? ~ ran ,, i s almost a constant, 
ig. 24 also shows the variation in the maid mum film thick- 


ness at X = C, with respect to p mj _ n . As 


p min increases , 




c 111 d f ilO 11C G G t/ 


ris es , reducing the co nd ens ation rate 


vm 


„ u v, 


decreasing the flow rate. Alternative!; 


‘tain 5 


would be known if 


ma: 


is 


the film profile and, hence -,p 
knovm. 

As is evident from Eq., (3.36), the film profile is 
a stroiir function of the length of the condensing surface, L. 
,?ig. 25 shows the condensate film profiles for a range of 1. 

It is observed that 6 fflax increases as 1 is increased, which 
is ounce ted, since the driving force must be greater for 
greater L. Also a larger length of the condensing surface 
L vives a .great or condensation rate and, hence, a larger 
flow rate. 

The film profile and the rate of condensation are 
also affected by the properties of the condensing vapour 
pnd b th0 difference between vapour saturation temperature 
and surface temperature. At. Jig. 26 presents the condensate 
-film profiles for various values of X, defined in Eq.(3.24). 
The value of X increases with the viscous forces, the thermal 
conductivity and the temperature difference. At, and 
doorcases when gravity forces or latent heat of the eoixdensin £ 
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^rrvDOur inc reason. And , hence, the value of S is large for 
a highly viscous condensing fluid. Pig. 26 shows -chat the 

- ^ •? -f'v.--o S cs with K and, thus, the condensation 

w x s re cl. uo w * 

flic local Huso olt number IIu^, at various locations 
has Icon presented in Pig. 27, for a range of ri mirL keeping 
X r . n d 1 as constant. As r) rain decreases, hu x at xhe end of 
the condensing surface increases, while the variation in its . 
value at low values of X is quite small. The value of mean 
tjueoclb nuaber, :!<V for • *«» ’■' hoa ^ain = °' 4 ’ ls fom113 t0 

bo within 2% of 0.82 Sh l/5 M reported by Ha and 
Loppert [21] but the value is appreciably larger for 

"*) « £. 0 * *'• * 

*min 

The effect of surface tension terms in Eq.(3.40), for 

, . . d coords on the values of Bond number, 

steady, laminar flow, aepenew 

„„ , „ o-n^r-i-od to be larger for smaller values 

Bo. The effects are cspec vea to ° 


0 x u 


30. Eq. (3.40) lias been solved by finite - difference 

r. Trnlues of ri . and the length of the 

techniques lor various values ox ii^ 

„ t mi no results show that the effect of 

condensing surface, L. The igsux^ 

.-r. -nnp-lio-ible, even for very small values 
surface tension is nc&iXb 1 ^’ 

of the Bond number. The change in maximum film thickness, 

5 , is within 0.1 percent for Bo > 0.0001, when u min 

equals 0.6 and 1 is 4 ems. The effect of surface tension 
diminishes, largely due to .assumed round fall at the end 
Of the condensing surface. Per the condensing surface 




j w } 

- f o 1 o 

cx>i N| in 





th 


+4 


tormina ted by a sharp edge, tile effect is expected to be 
cons icier : ole , but ouo to ties lack ox knowledge of proper 
bound ~ ry conditions, it could not be considered. 

(b) Transient , lam innr Film Condensation Flow: 

ic is clour from the Eqs, (3.22) and ( 3 . 27 ) , the 
transient effects on condensate film profile are controlled 

given by Eq. (3,27). The parameter R is 
btainod by Sparrow and Siegel [18] for 
transient laminar film condensation on a 
The values of R depend on the condensing 
length of the condensing surface, 1, and 
by, g and, as such, the transient 
r highly viscous condensing fluids, 
vity. Pig. 28 presents the condensate 
film profile, with negligible surface tension, for various time 
obtained for the initial condition that the film thickness is 
uniform at % = 0 and , after the flow starts, b min is a cons- 
tant. It is observed that the rate of increase in film 
thickness doorcases with time, and, hence, the local heat 
transfer coefficient h„_, also decreases. The effect of 

JL 

surface tension is expected to be the same, as discussed fox 1 
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the steady, laminar flow 
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-f ^ y. 
U 'w- J. 



( b ) 


minuted jy a sharp edge, the ef 
siderr ole , but duo to tho lack 
ndary conditions, it could not 


transient , Lara in 


Pi In 


^ O 4-i' 


d on 


foot is expected to bo 
of knowledge of proper 
be considered. 

cation Plow: 


-~- c xs clear from the Eqs, (3.2 2) and ( 3 . 27 ), the 


t ms i o r. t off 


by a parameter R, given by Ea. (3.27). Tho parameter R is 


tho same as 


that obtained by Sparrow and Siogcl [18] for 
their ana .lysis of transient laminar film condensation on a 
vortical surface . The values of R depend on the condensing 
vapour properties, length of tho condensing surface, 1, and 
acceleration duo to gravity, g and, as such, tho transient 
effects will be larger for highly viscous condensing fluids, 
or for a field of low gravity. Pig. 28 presents the condensate 
film profile, with negligible surface tension, for various timer 


obtained for the initial condition that the film thickness is 
uniform at x = 0 and, after the flow starts, p .. is a cons- 
tant. It is observed that the rate of increase in film 
thickness decreases with time, and, hence, the local heat 

transfer coefficient h , also decreases. The effect of 

X 

surface tension is expected to be the same, as discussed for 
the steady, laminar flow. 



CONCLUSI OilS 


CHAPTER IV 


The results of the present investigation are import. ant 
as much of this work is the first known treatment of many of 
the problems. for the transient film condensation on the 
upper surface of a horizontal plate, which is considered foi 
five different boundary conditions, it is emphasized that 
there are many possible physical mechanisms which could govern 
run-off. The analysis is done for two inviscid flow conditions, 
one, a specified fraction of the rate of condensation, 
occurring at any time, flowing out across the edges of the 
condensing surface at that instant, and, second, the out- 
flow dominated by gravity. The theoretical results are 
shown to depend on two dimensionless parameters, one for con- 
densation, p, except for the case of iso thermal surface, and 
the other for the flow, r or y. The situation of no-flow 
arises when the flow parameter is zero, and, hence, transient film 
condensation . with no run-off depends on only ono dimension- 
less parameter, p. This result would be applicable for short 
times when the flow has not started due to surface tension 
effects and for circumstances where a wall obstructs bhe 

condensate flow. 
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Lami na r f i In c o ncl one ati o n f 1 o w on a- h 0 ri 0 n lj a -*- 
iso thermal aurfa.ee is analysed in three stages, (a) transient 
film condensation without flow, (b) transient, l&minaj. flow 

and (c) steady, laminar flow. The analyses for cranoimt 

, , „ , , ..^face tension effects 

and s*ceao.y laminar flows consider the 

_ . , . , . „ „ eauations are 

also, fourth order non— linear duf ore nr o — 

, , . , „ ,, . , . , +ho transient and 

obtained for the xilm thickness m been 

_ „ ...^n<?ion. or large values 

the laminar eases, negligible surface 

„ ^ , ,, , ^^anr systems. The 

of Bond number reduce them to second oro-i 

steady, laminar condensation flow is shown depend upo 

parameter, K, which is a function of vapour properties 

, , , , the vapour and 

the difference between the temperatures ox 

m , , . _ f-iim profile for the 

the condensing surface. The condensate xj-j- 

above flow is found to bo a strong function of the 

_ with length of 

condition and the film thickness increase^ 

, - + W K. The transient 

the condensing surface, and the parameter 

, . , i ed by another 

response of the film thickness is controixou j 

, „ , , 0 vaxour properties, 

parameter R, which, also, depends on une 

j_v, ~ acceleration duo to 
length of the condensing surface and tne - 

_ , , TT the transient 

gravity. The analysis and results show on 

a pns ins fluids and 

effects arc large for highly viscous con 

-i that for oho 

for low gravity fields. It is also oboer 

-fall, the surface 

condensing surface terminated by a roun<-- 

„ vprV small values of 
tension effects are negligible oven for 


Bond number. 



For th transient film c one! e nsat ion , conoinea wibh 


conauction, on a horizontal 


the future work may be 


concentrated in the following areas ; theory say bo developed 


for the flow fraction pa raneter r as a function 01 ci'.i 


and 


e up crime ires need, to be performed on substantiating the 
theory d eveJ.oped . 

For the laminar film condensation flow on a horizon- 
tal isothermal surface, further work is needed, on several 
aspects . 

For transient film condensation, without flow, sufficient 
information must be obtained for the conditions jus u before 
the start of the flow, for various edge conditions . 

By experimentation or analysis, information should oe 
obtained for specifying the boundary condition a'c che edge 
of the condensing surface, terminated by a sharp edge, m 
the absence of which transient end steady state soluuions 
could not be obtained for the sharp edge ca.se. The efxeci.s 
of R and K need to be further examined for transient flow. 
The effect of surfaces tension, particular^ ~ or a s ^- ar P 
edged surfa.ee, needs to be studied further over a range 


of Bond number. 


Hussolt , 7. , 'Die Oberflachon Kondcnsation de 
Mass erdampf os' , Eoitschrift des Veroines deutcher 
Ingenioure, Vol. 60, 1516, pp. 541-546. 

Bromley, L.A. , 'Effect of Heat Capacity of Condensate 
in Condensing' , Ind. Engg* Chom., Vol. 44, 1552, 

o r\ r r no /Ca 

pp. 


ASME . 


Eons enow, 7.M., 'Heat transfer and. Temperature 
Distribution in laminar Film Condensation' , Trans 
Vol. 73, 1556, pp. 1645-1648. 

Sparrow, E.M. and J.L. Gregg, 'A Boundary-layer 
Treatment of Laminar Film Condensation', J. Heat Tr. , 

Vol. 81, 1559, pp. 13-18. 

Snarrou , E.H. and J.L. Gregg, 'Laminar Condensation 
Heat Transfer on a Horizontal Cylinder', J.Hcat Tr. , 

Vol. 81, 1955, PP- 251-256. 

Gebliart , B. , 'Heat Transfer', 2nd Edition, McGrow 
Hill Book Co., 1571, p- 443. 

Soli, J.C.T., E.M. Sparrow and J.P. Hartnett, ' Tho Two 
Phase Boundary Layer in Laminar Film Condensation , 

Int. J. Heat Mass Tr. , Vol. 2, 1961, pp. 6g-o2. 

Chen, M.M., 'An Analytical Study of Laminar Film 
Condensation: Part 1 - Flat Plates', J. Heat Tr. , 

Vol. 81, 1961, pp. 48-54. 

Poets, G. and R.G. Miles, 'Effects of Variable Physical 
Properties on Laminar Film Condensation of Saturated 
St cam on a Vertical Flat Plate' , Xnt. J. Heat Mass 
Tr., Vol. 10, 1967, pp.1677— lo9^. 

y anP - K.T., 'Laminar Film Condensation on a^Vertical 
Ho niso thermal Plate', J. Appl. Meeh., Vol. 33,1966, 
pp. 203-205. 

Sparrow, E.M. and S.H. Lin, 'Condensation Heat 
Transfer in the Presence of a Eon-condensable Ga- , 

J. Heat Tr, , Vol. 86, 19 o4, PP* 430—436. 

Minkowvcz, W.J* and E.M. Sparrow ,' Condensation Heat 
Transfer in the Presence of Non-condensables, Interfacral 
Resistance, Super Heating, Variable Properties, and 
Diffusion', Int, J. Heat Mass Tr. , Vol. 9, 1966, 
pp. 1125-1144* 



13. 


Sukhatmc 
luring ] 
J. Ho at 


•/ 

’ilr 


arxc; 


H T-; 


Rohscncw , f Heat 


; rano i or 


rn ^ 

* 9 


.,n Condon sat ion of a liquid Metal Van our ' , 

Vol, 83, IS 56, pp. 19-28. 


14. 


Mills 

of Ha’ 

1827. 


A.?, and R.A 

J . Heat 


Seban, 


' The 


1X1 G 


■’■i a qq J , *• Y* 
r:Cv wJ O 4_ X • 


C o nd one at io n C o o f f ic i ent 
• Vol. 10, 1967, pp . 1815- 


15. 

16 . 


17. 


1 


IS. 


20. 

21 . 


22 


23a. 


23b * 
24. 


ICoh, J.C.Y., 'film Condon sat io 
Boundary layer Plow', int. J. 
1S62, p P . 941-954. 


n in a Pore ad 
Ho at Hass Tr. 


Convection 

Vol. 5 , 


Sparrow , B.I*. , H.J. Hinkowycz and. 17. Saddy 
Convection Condensation in the Presence of 
sables and Int or facial Resistance' , Int. J 
Tr. , Vol. 10, 1967, pp. 1825-1045. 


1 Forced 
iTon-conden- 
Heat Mass 


Hinkowycz, H.J. and S.M. Sparrow, 'The Effect of 
Superheating on Condensation Heat Transfer in a 
Forced Convection Boundary Layer Plow,' Int. J.Heat 
Mass Tr. , Vol. 12, 1965, pp. 147-154. 


Sparrow, S.?I. and R. Siegel, 'Transient Film Condensation’ 
J. Anpl. Moch., Vol. 26, 'IS 59, PP. 120-121. 

G-ebhart, 3., 'Random Convection Under Conditions of 
Weightlessness', AIAA Journal, Vol. 1, 1963, pp. 380-333- 


Leppert, Of. and B. Fimno, 'Laminar Film Condensation on 
Surfaces Formal to Body or Inertial Forces', J. Hoat Tr. , 
Vol. 30, 1963, pp. 178-179. 

Himrao, B. and G-. leppert, 'laminar Film Condensation 
on a Finite Horizontal Surfo.ce', Heat Transfer 1970, 

Vol. VI, Elsevier Publishing Co., Amsterdam, 1970, 

Cs 2.2. 


Pfahl, R.C., Jr., 
Condensate Flow' , 
pp. 738r-793 . 


1 T ran s i en t Co nd on s a t io n Hit hou t 
Int, J. Heat Mass Tr. , Vol. 17, 1974, 


Rohs enow , 
Transfer, 


U„rl. and J.P. Hartnett, 'Hand Book of Hoat 
HcG-raw Hill Book Co., 1373, Chap. 12, p. 16. 


Rohsonow, U.M. and J.P. Hartnett, 'Hand Book of Heat 
Transfer, McGraw Hill Book Co., 1973, Chap. 12, p.8. 

Holland orf , J.C. , and T.Y. Chu, ’Transient Film 
Condensation on Finite Horizontal Surfaces with 
Run-Off, A SMI Paper 75-HT-16, 1376. 



Bankoff » S.G-. , 1 She Cento rtional Energy 
the Spreading of Large Drops 1 , J. Phys. 
1956,” pp. 952-355. 


Re quir cm e nt s i n 
Chem. , Vol. 60, 


Cars law , 
Solids , 

1373 , PP 


H.S. and J.C. 
Oxford Univers 
285 “283. 


Jaeger s 
ity Press 


'Conduction of Heat in 
, London, 2nd Edition, 


Goodman, 

Transient 

Transfer* 


JL 


. R. , 'Application of Integral Methods to 
Nonlinear Heat Transfer', Advances in Heat 
Vol. 1, Academic Press, New York, 1964, PP* 53-64 


Goodman, T.R. 
Heat Inputs ' , 


'The Heatinv of Slabs with Arbitrary 
J. Aerospace "Sci, , Vol. 26, 1959, pp. 187-188 


Hreith, P. . ! Principles of Heat Transfer', 2nd Edition, 
Int. Textbook Co., '1965, p. 128. 


Ingcrsoll, L.R. , O.J. 
Conduction, Indian Ed 
1969 , PP- 194-198. 


Sobol, and A.C. Ingersoll, ’ Heat 
, Oxford and IBII Publishing Co., 


Eckert, E.R.G. and R.M. Drake, Jr., 'Analysis of Heat 
and Mass Transfer* Int. Student Edition* McGraw 
Hogakusha, 1972, pp. 183-186. 


ICutateladse, S. 8. Fundamentals of Heat Transfer', 
Edward Arnold Ltd., London, 1963, pp. 300-307. 


Merte, II., Jr., 'Condensation Heat 
in Heat Transfer, Vol. 9, Academic 


Transfer', Advances 
Press, 1973, . p.186. 



A54885 


Pate Slip I; 4 8 S 


This book is to be returned on the 
date last stamped. 



CD 6.72.9 



H0£- F'L 



